Violation of Luttinger's Theorem in the Two-Dimensional t-J Model 
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We have calculated the high temperature series for the momentum distribution function rik of the 
2D t-J model to twelfth order in inverse temperature. By extrapolating the series to T — 0.2J we 
searched for a Fermi surface of the 2D t-J model. We find that three criteria used for estimating 
the location of a Fermi surface violate Luttinger's Theorem, implying the t-J model does not have 
an adiabatic connection to a non-interacting model. 
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Models for two-dimensional strongly correlated elec- 
trons play a central role in attempts to understand high 
temperature superconductors . However, the 2D mod- 
els themselves are at present poorly understood. One of 
the main points of interest in studies of 2D strongly corre- 
lated electrons is how similar the 2D models are to Fermi 
liquid theory, the standard model for conventional metals 
Many-body calculations for conventional metals are 
generally perturbative, assuming an adiabatic relation to 
a non-interacting model, with low energy excitations de- 
scribable by quasiparticles. 

By summing a perturbative expansion to all orders 
Luttinger |3| was able to show that a sharp Fermi surface 
can exist for interacting electrons. He defined the Fermi 
surface to be the locus of points in k-space where the 
renormalized single particle energy is equal to the zero 
temperature chemical potential = /U. This requires 
the imaginary part of the retarded self-energy to vanish 
on the Fermi surface. Luttinger was able to show Q] 
ImE ret (w) oc (uj — /i) 2 , satisfying this requirement. An 
immediate consequence of this perturbative calculation 
is that the volumes (areas in 2D) enclosed by the inter- 
acting and non-interacting Fermi surfaces are the same, 
a statement generally known as Luttinger's Theorem. 

Using high temperature series we investigated the mo- 
mentum distribution function for the 2D t-J model on 
a square lattice, with the Hamiltonian for the t-J model 
given by 



H 



C j<7 Ci, 



P 
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where the sums are over pairs of nearest neighbor 
sites and the projection operators P eliminate from the 
Hilbert space states with doubly occupied sites. We cal- 
culated the high temperature series for the momentum 
distribution function to twelfth order in inverse temper- 
ature (3 = 1/ksT, extending a previous eighth order cal- 
culation by Singh and Glenister || . The definition of the 
single spin momentum distribution function is 



n k 



(2) 



with n r = {c) 0(y c r(y ). For the calculations reported here 
we fix J/t = 0.4 and the electron density n — 0.8. A 
wider range of parameters will be explored in a future 
publication. 

To reach low temperatures we need to analytically con- 
tinue the series for nk- A standard way to do this is to 
use Pade approximants. However, for nk the straightfor- 
ward application of Pades does not work very well. One 
way to improve the convergence of Pade approximants is 
to change the functional form before calculating Pades 
||. Exactly which change to make is difficult to know 
for unknown functions. One way to proceed is to use 
the function itself as a scaling function. To do this we 
first form the high temperature series for the ratio of two 
values of nk with the two k-points closely spaced and all 
other parameters the same. This is an exact calculation 
using the exact coefficients we have for n^. The series 
for the ratio is then analytically continued using Pade 
approximants down to T = 0.2J (for J = 1500 K this is 
300 K). At each temperature all of the ratios are refer- 
enced back to the zone center with an absolute scale for 
rik set by enforcing the sum rule X)k nk = n /2- Using 
this technique we collected data for rik at 1326 k-points 
in the irreducible wedge of the square Brillouin zone for 
five temperatures T/J = 0.2, 0.4, 0.6, 0.8 and 1.0. The 
k-points have a uniform spacing of 7r/50 for k x and k y . 
The accuracy of the calculation varies with k, but for 
T = 0.2J is approximately 5%, for T — 0.4J approxi- 
mately 3% and for higher temperatures 1% or less. 

In analyzing our results our main goals were to search 
for a possible Fermi surface of the 2D t-J model and 
compare its area to the area of the tight-binding model 
Fermi surface to check Luttinger's Theorem. Since we 
cannot reach T = 0, where the Fermi surface is defined 
by a sharp discontinuity in rik, we consider three criteria 
which can be calculated for non-zero temperature and 
might be expected to smoothly approach the T = Fermi 
surface as T —> 0. The curves in k-space we investigated 
are defined by i) nk = 1/2, ii) dn^/dT — and iii) |Vk«k| 
maximal. 

A comparison of the locus in k-space where rtk = 1/2 



1 



for the t-J model at T — 0.2J to the tight-binding 
model Fermi surface is shown in Fig. 1. The area en- 
closed by the t-J model curve is smaller than the area 
of the tight-binding model Fermi surface. This differ- 
ence is not due to non-zero temperature. Fig. 2 shows 
the temperature dependence of the t-J model nk for 
ki? = (0.42667T, 0.42667r), the Fermi momentum of the 
tight-binding model on the zone diagonal. Clearly nk for 
the t-J model with this momentum remains below 1 /2 for 
all temperatures. The insert in Fig. 2 shows the tempera- 
ture dependence of the chemical potential /z(T) for the t- 
J model with the same parameters as for n^. The chemi- 
cal potential shows little variation for T ^ J, appropriate 
for a degenerate Fermi system, with fi k 1.66i = 4. 15 J. 
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FIG. 1. Comparison of the tight-binding Fermi surface 
(dashed line) to the locus of k-points for the 2D t-J model 
where = 1/2 at T — 0.2 J (solid line), both at an electron 
density of n = 0.8. The area enclosed by nk = 1/2 for the t-J 
model does not satisfy Luttinger's Theorem. 

This comparison alone is not sufficient to claim a vio- 
lation of Luttinger's Theorem. In the absence of particle- 
hole symmetry nk need not equal 1/2 on the Fermi sur- 
face raj. The 7ik = 1/2 criterion for the Fermi surface has 
been widely applied in the past || and is the simplest cri- 
terion to check before doing more detailed calculations. 
Our k- resolution allows us to distinguish nk = 1/2 for 
the t-J model from the Fermi surface of the tight-binding 
model with approximately three data points separating 
the two curves. Previous calculations on small clusters 
did not have sufficient k-resolution to make this dis- 
tinction. 

An improved criterion for locating a Fermi surface has 
been proposed by Randeria et al. Jj). They proposed 
that the Fermi surface be identified with the locus of k- 
points where the temperature derivative of nk is station- 
ary, drik/dT = 0. While the high temperature series for 
dn^/clT can be calculated directly from the series for nk, 
the resulting series is too short to be extrapolated to low 
temperatures. Alternatively, we consider the finite dif- 
ference approximation Ank/AT where AT = 0.2 J and 
Arik is found by directly subtracting the two rik's for 
each k-point. Fig. 3 shows that this works very well for 
the tight-binding model at T — 0.3J. From our data 



we find Ank/AT centered on the average temperatures 
T/J = 0.3, 0.5, 0.7 and 0.9. The results are shown in 
Fig. 4. 
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FIG. 2. Temperature dependence of the t-J model mo- 
mentum distribution function rik at the tight-binding Fermi 
wave vector kF = (0.4266-7T, 0.42667r) showing that on the 
tight-binding model Fermi surface the t-J model rik < 1/2 
for all temperatures. Insert: Temperature dependence of the 
t-J model chemical potential. For T ^ 0.5J the error bars 
are the width of the line. 

The results shown in Fig. 4 have three surprising fea- 
tures i) the area enclosed by the curve where Arik/AT = 
is larger than the area enclosed by the tight-binding 
Fermi surface, ii) k-states with decreased occupancy as 
the temperature increases occur across a broad region in 
the center of the Brillouin zone and iii) at low temper- 
atures the distribution of k-states with increased occu- 
pancy as the temperature increases are strongly peaked 
on the zone diagonal. Features ii) and iii) indicate the lo- 
cation of the low energy excitations in the Brillouin zone. 
Unlike the tight-binding model, where the low energy ex- 
citations are confined to a narrow range of k-points cen- 
tered on the Fermi surface, the low energy excitations for 
the t-J model are spread throughout the Brillouin zone. 
The extremal values of Ank/AT for the t-J model are 
~ 10 times smaller than for the tight-binding model, im- 
plying a relatively small density of states at any fixed mo- 
mentum for the t-J model. For dn^/dT = to give the 
correct Fermi surface particle-hole symmetry is required 
for low energies R|. This may not be valid for the 2D t-J 
model. Thus the locus of k-points where dn^/dT = 
may not give the true Fermi surface. 

The final criterion we considered for locating a Fermi 
surface for the 2D t-J model is to follow the locus of 
k-points where |Vknk| is maximal [@J^]. This criterion 
only depends on nk having a sharp discontinuity at the 
Fermi surface for T = 0. For T > the discontinuity 
will be smeared out, but we still expect |Vknk| to be 
large near the Fermi surface since we have a degenerate 
Fermi system. As shown in Fig. 3 this criterion gives the 
proper Fermi surface for the tight-binding model. For 
the t-J model we calculate |Vknk| numerically from nk, 
with the results shown in Fig. 5. 
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FIG. 3. Full Brillouin zone plots for the tight-binding 
model. Left: An k /AT for f = 0.3J = 0.12t, Right: |V k n k 
for T = 0.2J = 0.08t. For both plots blue is zero. In 
the left plot the orange, yellow and green regions are neg- 
ative with a minimum value of An k /AT = —0.9 J -1 and the 
red and violet regions are positive with a maximum value of 
An k /AT = 0.6J~ . In the right plot red and violet regions 
are positive with a maximum value of |V k n k | = 6.8a, where 
a is the lattice spacing, here set to one. Also, in both plots 
the solid black curve is the Fermi surface of the tight-binding 
model, showing excellent agreement with both approximate 
techniques for locating the Fermi surface. 




tight-binding model, the maximum value of |Vk^k| for 
the t-J model at T = 0.2J is ~ 10 times smaller. The 
shape of the area enclosed by the dotted line in Fig. 
5 is similar to that for the tight-binding model with a 
next-nearest neighbor hopping term, but this modifica- 
tion to the tight-binding Hamiltonian can only change 
the shape of the Fermi surface, not its area. Thus our 
results for the t-J model cannot be modeled as a band of 
non-interacting electrons with a modified band structure. 




FIG. 5. Full Brillouin zone plots of |V k n k | for the t-J 
model. All of the colors are values of |V k n k | > 0, with orange 
in the zone center, around k = (n, n) and around k = (n, 0) 
being zero. Red is the maximum value V k n k | = 0.43a, where 
a is the lattice spacing, here set to one. The plots are for 
different temperatures. Upper left: T = 0.8J, Upper right: 
T = 0.6J, Lower left: T = 0.4J and Lower right: T = 0.2J. 
In the lower right plot the solid line is the tight-binding Fermi 
surface, the dashed line is the locus of k-points from Fig. 
4 where An k /AT = and the dotted line is the locus of 
k-points where |V k n k | is maximal (only one of four branches 
is shown). 



FIG. 4. Full Brillouin zone plots of An k /AT for the t-J 
model. The color scale is the same for all four plots. Or- 
ange, yellow and green are negative, with a minimum value 
of Ank/AT = -0.08 J -1 . Blue, violet and red are positive, 
with a maximum value of An k /AT = 0.15J . The plots 
correspond to different temperatures. Upper left: T = 0.9J, 
Upper right: T = 0.7J, Lower left: T = 0.5J and Lower right: 
T = 0.3J. In the lower right plot the solid black curve is the 
tight-binding model Fermi surface and the dashed curve is the 
locus of k-points where An k /AT = for T = 0.3 J. 

At the lowest temperature shown in Fig. 5, T = 0.2J, 
our results have two main features i) the area enclosed 
by following a continuous locus of k-points along a ridge 
where | Vk^kl is maximal (as indicated by the dotted line 
in Fig. 5) is larger than the area enclosed by the tight- 
binding model Fermi surface and ii) |Vk^k| is strongly 
peaked on the zone diagonal. However, compared to the 



The main result of our calculation is that the 2D t- 
J model violates Luttinger's Theorem. This means the 
ground state of the 2D t-J model is not adiabatically re- 
lated to a non-interacting model. Also, the distribution 
of low energy excitations revealed by Ank/AT suggests 
that quasiparticles cannot describe all of the low energy 
degrees of freedom of the 2D t-J model. However, our re- 
sults are not sufficient to uniquely determine the ground 
state of the 2D t-J model. 

The simplest way to not have an adiabatic connection 
to non-interacting electrons is for the 2D t-J model to 
have an ordered ground state with a different symme- 
try than non- interacting electrons. The low temperature 
growth of peaks in Ank/AT and |Vk»ik| suggests some 
kind of order is developing in the t-J model. The entropy 
of the t-J model [[To) also starts to decrease at T ~ J. At 
present it is not possible to determine the precise nature 
of this order. We know from the spin correlation function 
pT| and the antiferromagnetic correlation length |12j that 
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we do not have long range spin order. Concomitantly, we 
do not observe hole pockets in rik Jl^|. For the t-J model 
the charge fluctuations are suppressed, with no indica- 
tion of short wavelength structure in the charge density 
jyl]. The location of peaks in Ant/AT and | Vk»ik| along 
the zone diagonal is consistent with the location of gap 
nodes for d x 2_ y 2 superconducting fluctuations. If this 
is correct, superconducting fluctuations at low tempera- 
tures are developing in the 2D t-J model from a higher 
temperature state which cannot be described as a Fermi 
liquid. For T ~ J the t-J can be modeled by assum- 
ing spin-charge separation |l4|,[l5|. However, there is no 
clear reason for peaks to appear in An^/AT or Vk«k| 
for a spin-charge separated state. 

Other theoretical approaches have found results simi- 
lar to our data. In particular SU(2) gauge theory calcu- 
lations |16| and phenomenological models of preformed 
pairs p!7| or of a d-wave ground state pl| find Fermi 
arcs centered on the zone diagonal and the Fermi surface 
gapped or destroyed by strong scattering near k = (n, 0). 
These theories are motivated by the pseudogap observed 
in ARPES @. The largest value of |V k rt k | and the 
sharpest region for Ank/AT = form arcs centered on 
the zone diagonal with the same outward curvature as the 
ARPES data. However, the maximum value of |V k ^k| is 
substatially smaller than what might be expected for a 
true discontinuity unless there is an energy scale in the 
t-J model smaller than 0.2 J = 300K. Our data also are 
similar to the "cold spots" proposed by Ioffe and Mil- 
lis [|o| and quasiparticle decay as discussed by Laughlin 
pl| , plus calculations for t-J ladders that show gapped 
and non-gapped features in |p2f . 

ARPES experiments do not directly measure n^, 
though attempts have been made to integrate the 
ARPES data over frequency, giving a value proportional 
to nk, but with an unknown scale factor Q[. We can make 
a qualitative comparison of our results to features of the 
ARPES data |2j| i) |Vk«k| maximal gives a locus in k- 
space with the same shape as the Fermi surface found in 
ARPES experiments, ii) without Luttinger's Theorem, 
the peaks observed by ARPES in the normal state JtJ do 
not have to be quasiparticles and iii) we see a large re- 
gion with low energy excitations, similar to the flat back- 
ground extending up to the chemical potential in ARPES 
data 

In conclusion, we find that the t-J model violates Lut- 
tinger's Theorem. Using |Vknk| as the least biased way 
to search for a Fermi surface we find the area enclosed by 
the curve where |Vk?ik| is maximal to be larger than the 
area enclosed by the Fermi surface of the tight-binding 
model. We also found structure in the low energy excita- 
tions of the 2D t-J model which might be due to d x i_ y i 
superconducting fluctuations. 
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